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A study of emission temperatures extracted from excited state populations 
and of freeze out radii from light particle intensity interferometry is presented. 
Three high resolution Ai?~£'-Hodoscopes with a total of 216 detectors are 
combined with the ALADiN setup in order to study Au+Au collisions at 1 
A-GeV. In contrast to measurements with the isotope thermometer the ex- 
tracted apparent temperatures do not vary with impact parameter thus with 
excitation energy. From the extracted radii a freeze out density was determined 
which decreases from 0.2po for the most peripheral to less than O.lpo for the 
most central collisions. A density-dependent feeding correction is applied to 
the different temperature measurements. 
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1 Introduction 



In a recent publication |I| we presented the determination of emission tempera- 
tures and excitation energies in Au+Au collisions in order to explore the liquid-gas 
phase transition of nuclear matter. One open question was the influence of feed- 
ing on the isotope thermometer TneLi- In a follow up experiment we compared the 
isotope thermometer with temperatures derived from the relative yields of excited 
particle-unbound states. The idea was to calibrate the isotope thermometer with a 
temperature determined from the excited states rather than deriving this calibration 
from theoretical calculations. 

The following excited state thermometers were analyzed: 

1. ^Li which decays into (i-'^He {E* = 16.66 MeV) and into p-a in the ground 
state. 

2. ^Li where the apparent temperature was determined by the ratio of the third 
{E* = 4.31 MeV) and fifth {E* = 5.65 MeV) excited states to the first {E* = 
2.168 MeV) excited state which all decay into d-a. 

3. ®Be where the ratio of the state at E* = 17.64 MeV (decays into j9-^Li) to 
the state at E* = 3.04 MeV (decays into a-a) is used for the temperature 
determination. 

The paper is organized as follows: 

In section ^ details of the experimental setup are given. The extracted temperatures 
for various event selection criteria are presented in section ^ Extracted radii and 
freeze out densities from p-p and d-a correlations are presented in section ^. In 
section ^ the results are compared to feeding calculations. 

2 Experimental setup and data analysis 

The experiment was performed at the SIS accelerator at GSI, Darmstadt. The ex- 
perimental setup ist shown in Fig. ^ A 25 mg/cm^ Au target was irradiated by a 
1 A-GeV Au beam at a beam intensity of 10^ particles per second. The forward 
hodoscope ZDO consisted of 36 phoswich detectors with a coverage of 6.5< 9 <21.5'' 
and 0< < 360*^. Three small-angle high-resolution AE — i?-hodoscopes with a 
total of 216 elements were placed at laboratory angles 6 = 100 — 150°. Each of the 
216 elements consisted of a silicon detector of 300 yum thickness and a cesium iodide 
detector of 6 cm length for 160 detectors and 10 cm length for 56 detectors, read out 
by photodiodes. Isotopes were resolved for Z = 1 — 4. The energy calibration was 
performed according to ^ and an energy resolution of about 1% was achieved. From 
the cross calibration between detectors we got an additional smearing of 1% resulting 
in an overall uncertainty of ~ 2%. The superposition of the AE-E distributions of 
all detectors is shown in Fig. ^ The separation of the Hydrogen isotopes (at Z = 1), 
the Helium isotopes (at Z = 2) and the Lithium isotopes (at Z = 3) is clearly visible. 
The trigger required two valid hits in the ZDO and two valid hits in the combination 
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Figure 1: Experimental setup 



of the three hodoscopes. For event characterisation we used the ALADiN Time-of- 
fiight wall, which is described in detail elsewhere [Q. The variable Zi,ound in the TOF 
array which is the summed charge in forward angles excluding the Z = 1 particles 
serves well as a criteria of the centrality of a collision Analogous event selection 
criteria as in |1|, Q were used. By the measurement of the decaying projectile specta- 
tor fragment charges with Z > 2i and the neutrons (0) an excitation energy E* and 
a mass of the prefragment Aq could be assigned to every Zbound bin [Q. The mass of 
the prefragment Aq was used in the present analysis to extract the freeze out density 
from the radii derived by p-p and d-a correlations. 

In Fig. ^ energy spectra for light particles together with moving source fits assuming 
three sources in a symmetric system are presented. The moving source fits gives the 
following results: 

1. The main contribution of the spectra can be described by a single source, the 
target source moving with less than 3% c in the laboratory frame. 

2. The inverse slopes of the spectra decrease with increasing mass of the particle 
species even for the most central collisions. 

3. The inverse slopes increase continuously with increasing centrality. This rise is 
partly due to the increase of bounce with decreasing impact parameter partly to 
the increasing excitation energy with decreasing impact parameter. Its difference to 
temperatures extracted from isotope yields or excited states can be partly assigned 
to the effect of the Fermi motion ||^ and to a significant part of pre-breakup emission 
of light particles . 
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In order to construct the correlation function for a specific particle pair one combines 
all particles of one species with all particles of the other species within one event and 
sorts the pairs according to their relative momenta. The coincidence spectrum shows 
the resonances of interest on top of a broad background distribution which is de- 
termined by phase space and acceptance of the hodoscopes. This background was 
determined with the event mixing technique combining particles from different 
events. The correlation function is constructed by dividing the coincidence yield by 
the background yield and normalizing it to 1 for large relative momenta as shown in 
Fig. ^ for the d-a correlation function. For the determination of a temperature one 
is interested in the absolute yields of the different states. First, the nonresonant part 
of the correlation function has to be subtracted. This part comes from Coulomb re- 
pulsion and quantum mechanical effects and was determined by scaling non-resonant 
correlation functions like t-t or t-^He. After subtracting this non resonant part 
of the correlation function and multiplying the resulting spectrum with the integral 
coincidence yield a resonance spectrum is derived. This is shown in the lower right 
panel of Fig. ^ The spectrum is compared to a Monte Carlo calculation for a certain 
emission temperature (here 5 MeV), which includes the first five excited states of ^Li, 
the detection efficiency and resolution of the detector. Details of the Monte Carlo 
calculation are presented in Ref. [§]. 



3 Temperature extraction from excited states 

Details of the level schemes for the analyzed thermometers are presented in table |^ 
For the ^Li thermometer the ratio of the first excited state at E* =16.66 MeV to 



Figure 2: AE vs E (left), PID vs energy in MeV (middle) and PID spectrum (right) 
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Energy spectra and moving source fits 
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Figure 3: Energy spectra at (6) = 35° and 2 source fit for 1 A-GeV 
dashed lines = target source, dot dashed lines = midrapidity source 
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Table 1: Level schemes of analyzed particle-unstable states (from 0) 

the ground state was used. The wide ground state was fitted with the R-matrix 
formalism with values taken from 0. The two correlation functions together with 
a minimum and a maximum background are shown in the first row of Fig. ^. The 
well known peak of ^B— >• p+^Be^.^ ^ p + {a + a) just above threshold could not 
be resolved in the experimental correlation function. The apparent temperature for 
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Figure 4: Coincidence yield in the d-a correlation function (upper left), uncoincident 
yield from the event mixing technique (upper right), ratio of coincident and unco- 
incident yield (symbols) and assumed nonresonant part (line, lower left), extracted 
yield (symbols) compared with simulated yield (line, lower right) 



^Li and the other excited state thermometers is determined by fitting the yield ratios 
of simulated resonance spectra for a range of input temperatures to the experimental 
yield ratios and using to the exponential temperature formula 

^Be decays in the ground and the first excited state into a-a, the fourth to sixth ex- 
cited states decay into p+'^Li (Fig. |^ second row). The ground state (g = 18 MeV/c) 
is too close to the detection threshhold of g = 15 MeV/c, so that we used the ra- 
tio of the fourth excited state at E* = 17.64 MeV and the first excited state at 
E* = 3.04 MeV for the temperature determination. The disadvantage of this ther- 
mometer is a rather small yield of ^Be resulting in large statistical uncertainties. 
Thus the apparent ^Be-temperature was determined without any centrality selection 
criterium. 

For the ^Li thermometer we have chosen the ratio of the 4.31 MeV and 5.65 MeV to 
the 2.17 MeV state. All three states decay into d-a (Fig. ^ lower right panel). Both 
resonances of ^Li were described with a Breit-Wigner function (see Fig. Because 
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1GeV excited state correlation functions 
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Figure 5: First row: Ground state (p-a, left) and first excited state (d-^He, right) 
of^Li 

Second row: Ground and first excited state (a-a, left), fourth to sixth excited state 
(p-Li, right) of^Be 

Third row: First, third and fifth excited state of (d-a, right) of^Li 



of the small energy difference strong feeding distortions of the Li population are 
expected. Additionally a small change in the the background correlation function 
results in a big change of the relative yields due to the large values of the background 
yields resulting in large systematic errors. No centrality selection criterium was ap- 
plied on the ^Li thermometer. 

The ^Li thermometer yields a rather constant apparent temperature of T 5 MeV 
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The apparent temperature values are in 

a- 



10, 11, 12, m 



independent of impact parameter (Fig 
the same range as reported elsewhere 

On the other hand, the extracted apparent temperature does not show a systematic 
rise as does the temperature Tuehi, which was derived from the relative isotope ratios 
of Helium and Lithium (open circles in Fig. 

For ^Li the dependence on the kinetic energy of the particle pair was also investigated 
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Figure 6: Apparent temperatures from excited states for different thermometers vs 
Abound and comparison with isotope thermometer TneLi from the target spectator 
without feeding correction — The error bars represent the statistical uncertainty, 
the systematic uncertainty is indicated by the brackets 



but no systematic trend could be established. 



4 Freeze out density from p—p and d—a correlations 

In the first row of Fig. U\ the p-p correlation function for four different cuts in 



Zbound together with a simulation from the Koonin-Pratt-formalism for Tq = 
8.0, 8.5, 9.0, 9.5 fm and r = fm/c is shown. 

The peak height of d-a and p-a correlations should also be sensitive on the freeze 
out radius [O, O- The dependence of those on Zbound is shown in the second and 



third row of Fig. ^ For the d-a correlations the integral of the first and second 
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Centrality dependence 




2.5 



Q 



2 



■A 



-A- 



4 



u- 

a- ■ 



o- - 



d-a 



rA- 



ill 




p-a 



50 100 50 



100 50 



0.3 r 



100 50 100 

q [MeV/c] 



0.2 



0.1 



r 



r 



PdH 



10 



30 



50 



70 



'bound 



Figure 7: First row: correlation function for four different cuts in Zbound (see 
X axis on fourtli row) and comparison to simulated correlation functions for Tq — 
8.0 — 9.5 fm respectively- 
Second row: d-a correlation function for tie same cuts as p-p in Zbound 
Third row: p-a correlation function for the same cuts as p-p in Zbound 
Fourth row: Mean density from extracted radii and prefragment mass Aq vs Zbound 
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maximum was compared to the integral of theoretical d-a correlation functions 
|11|, IIB, O in order to get absolute numbers for the d-a freeze out radius. The hard 



sphere radii vary from 10 fm for Zbound< 20 to 11 fm for Zbound> 60. 
Taking into account the increasing mass of the target spectator Aq with increasing 
impact parameter we determined a mean freeze out density. This mean freeze out 
density varies from 0.2po (0.17po) in peripheral to O.lpo (O.OTpo) in central collisions 
for the p-p (d-a) correlations (fourth row of Fig. |^. These densities are smaller 
than the standard density of 0.3po which is usually used in theoretical calculations. 
A freeze out density which decreases with decreasing impact parameter was already 
proposed by |l^ from INC calculations. 

In Fig. H p-p, d-a and p-a correlations are shown for for three different ranges of 
summed energy. The peak heights of the correlation functions and as a consequence 
the extracted radii vary strongly with the summed energy of the particle pairs (Fig. 
Faster particles seem to be emitted from a smaller source. This may be interpreted 
as a time ordered emission during expansion. 



5 Comparison to Feeding calculations 

For the analysis of feeding distortions we used the QSM code |T^. The Quantum 
Statistical Model is based on the assumption of thermal and chemical equilibrium. 
This model incorporates infinite system size. The ground and excited states of nuclear 
states up to Z = 13 are included according to published values 0. No continuous 
states are included in the QSM. The excluded volume of a fragment is the sum of 
the nucleon volumes (^H and ^He have the same excluded volume). The states are 
populated according to one global temperature T, density p and N/Z ratio. The 
initial particle ratios are modified due to the sequential decay of particle unbound 
states (feeding) and one has to distinguish between an initial and a final distribu- 
tion. The final ratio was used to calculate yields of ground and excited particle 
unbound states. By applying the exponential temperature formula (equation [l|) we 
derived a relationship between the initial temperature and the apparent ^Li, *Be and 
^Li temperature respectively. The comparison of the experimental yield with the 
QSM calculation gives the following results: 

1. The constant correction factor for the TneLi thermometer of 1.2 used by |l| 
was determined assuming a constant freeze out density p = po/3. However, the 
extracted freeze out radii from d-a correlations measured in the same experiment 
can best be described with a density varying from 0.17po for the most peripheral 
to 0.07po for the most central collisions. Because of the density dependence of the 
QSM correction the correction factor for the isotope thermometer becomes centrality 
dependent ranging from 1.2 for the most peripheral to almost 1 for the most central 
collisions. This reduces the observed effect of a steep rise in temperature for the most 
central collisions. The TncLi thermometer before and after QSM correction is shown 
in Fig. |. 



10 



Energy dependence 
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Figure 8: First row: correlation function for three different cuts in summed 
energy (see liorizontal error bars in bottom panel) and comparison to simulated 
correlation functions for tq = 9.0 fm, tq — 8.5 fm and Tq = 7.5 fm respectively 
Second row: d-a correlation function for tlie same cuts as p-p in summed energy 
Third row: p-a correlation function for tlie same cuts as p-p in summed energy 
Fourth row: Rise of the integrated maximum in the correlation function with the 
summed energy for p^, d-a and p-a correlations 
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Figure 9: TneLi before and after density dependent feeding correction according to 
QSM and comparison to apparent ^Li temperature 



2. According to QSM there is no feeding correction factor needed for the extraction 
of the ^Li temperature. Thus within QSM and the assumption of a global freeze out 
for all particle species only part of the discrepancy between excited state and isotope 
temperatures can be explained. 

3. Additional feeding calculations for the ^Li thermometer were performed analo- 

The correction for the ^Li thermometer from this code 



gous to |gO| (see also |£T 
is less than 1 MeV. 



6 Summary 

The extraction of apparent temperatures from of excited state population ratios for 
the system Au+Au at 1 A-GeV was presented and temperature values ofT^b MeV 
are derived. In contrast to the isotope thermometer presented in ^ the apparent 
temperature stays constant with decreasing impact parameter thus increasing excita- 
tion energy. The freeze out densities extracted by the means of d-a correlations and 
used as input for QSM calculations seem to favour a decreasing freeze out density 
with increasing excitation energy. After feeding correction one gets a rather slow rise 
in emission temperature starting at T = 5 MeV for the most peripheral and reaching 
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around T = 9 MeV for the most central collisions with the isotope thermometer while 
the excited state thermometer saturates around T ?s 6 MeV emission temperature. 



This discrepancy is even more pronounced in Ref. p2[ and explained by different 



freeze out times for isotope and excited state thermometers. 

References 

[I] J. Pochodzalla et al., Phys. Rev. Lett. 75 (1995) 1040. 
[2] C.J.W. Twenhofel, NIM B51 (1990) 58. 

[3] A. Schiittauf et al, Nucl. Phys. A607 (1996) 457. 

[4] C. Grofi, PhD thesis, J.W.Goethe-Universitat Frankfurt, 1997 (unpub- 
lished) . 

[5] W. Bauer, Phys. Rev. C51 (1995) 803. 

[6] W. Trautmann, contribution to this Proceedings 

[7] G.I. Kopylov, Phys. Lett. B50 (1974) 472. 

[8] V. Serfiing, PhD thesis, J.W.Goethe-Universitat Frankfurt, 1997 (unpub- 
hshed) . 

[9] F. Ajzenberg-Selove, 

Nucl. Phys. A413 (1984) 28. 
Nucl. Phys. A490 (1988) 1. 

[10] J. Pochodzalla et al., 

Phys. Rev. Lett. 55 (1985) 177. 
Phys. Lett. 161 B (1985) 275. 

[II] J. Pochodzalla et al., Phys. Rev. C35 (1987) 1695. 

[12] G.J. Kunde et al., Phys. Lett. B272 (1991) 202. 

[13] T.K. Nayak et al, Phys. Rev. C45 (1992) 132. 

[14] C. Schwarz et al., Phys. Rev. C48 (1993) 676. 

[15] S.E. Koonin, Phys. Lett. 70B (1977) 43. 

S. Pratt and M.B. Tsang, Phys. Rev. C36 (1987) 390. 

[16] D.H. Boal and J.C. Shillcock, Phys. Rev. C33 (1986) 549. 

[17] D.H. Boal, C.K. Gelbke an B.K. Jennings, Rev. of Mod. Phys.62 (1990) 553. 

[18] P. G. Warren et al., Nucl-Ex-96 10008 (1996). 

13 



[19] D. Hahn and H. Stocker, Nucl. Phys. A476 (1988) 718. 

[20] X. Hongfei, W.G. Lynch, M.B. Tsang, W.A. Priedmann, MSUCL-1040 
(1996). 

[21] B. Tsang, contribution to this Proceedings 
[22] C. Schwarz, contribution to this Proceedings 



14 
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